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We propose sub-harmonic resonant optical excitation with femtosecond lasers as a new method for 
the characterization of phononic and nanomechanical systems in the gigahertz to terahertz frequency 
range. This method is applied for the investigation of confined acoustic modes in a free-standing 
semiconductor membrane. By tuning the repetition rate of a femtosecond laser through a sub- 
harmonic of a mechanical resonance we amplify the mechanical amplitude, directly measure the 
linewidth with megahertz resolution, infer the lifetime of the coherently excited vibrational states, 
accurately determine the system's quality factor, and determine the amplitude of the mechanical 
motion with femtometer resolution. 



In recent years, nanophononic and nanomechanical 
systems have emerged as intriguing subjects for studying 
mechanics, heat transfer and opto-mechanical coupling 
on a nanometer scale [ll-Q. From a fundamental point 
of view, they provide a route to study mechanical 
excitations and their interactions with other elementary 
excitations P, ^From an applied perspective they 
have opened a pathway for high sensitivity sensors in 
the zeptogram mass range and in the attonewton force 
range [a, ISj ■ In established experimental methods these 
systems are driven electrically, magnetically, thermoelas- 
tically f^, via radiation pressure from continuous wave 
lasers |8i], or via other optical non-radiation-pressure- 
based schemes ■ The frequencies of typical systems 
investigated so far are in the megahertz to gigahertz 
frequency range [3] ■ The investigation of higher frequen- 
cies is strongly restricted by the driving and detection 
methods. Here, we report a new method for the investi- 
gation of a vibrational system by sub-harmonic resonant 
excitation with a high-repetition rate femtosecond laser. 
This excitation scheme can be regarded as tuning the 
separation of modes of the frequency comb of a femtosec- 
ond laser ll| to a commensurable of the frequency of 
the phononic system. By sweeping the comb spacing of 
the femtosecond laser, resonant impulsive excitation of 
the mechanical oscillator can be achieved, which allows 
the determination of its quality factor in the gigahertz 
to terahertz frequency range with femtometer sensi- 
tivity for the mechanical amplitude. We demonstrate 
the amplification of the fundamental cigcnmode of a 
free-standing silicon membrane at 19 GHz by a factor of 
20 compared to the off-resonant case and determine its 
quality factor. 



The dynamical properties of the free-standing sil- 
icon membranes were investigated by performing fs 
resolution pump-probe experiments using the recently 



12, m 



developed high-speed asynchronous optical sampling 
(ASOPS) method 



This method is based on 
two asynchronously linked femtosecond Ti: sapphire ring 
lasers of repetition rate ^ 1 GHz. One laser provides 
the pump beam and the second laser the probe beam. 
In this technique the time delay between pump-/probe- 
pulse pairs of the two pulse trains is realized through 
an actively stabilized 10 kHz repetition-rate-offset A/^ 
between the two lasers. This allows the stroboscopic 
stretching of the time between two consecutive pump 
pulses (Ins) to 100 /iS with an increment in time delay 
between pump and probe pulses of 10 fs. The time re- 
solved reflectivity change of the investigated system is 
recorded with a detector of ^ 100 MHz bandwidth. The 
absence of mechanical moving parts in this technique in 
combination with the extremely high scanning rate A/^ 
allows obtaining a signal-to-noise ratio above 10^ in a 
characteristic acquisition time of a few seconds [13, [3] ■ 
The limit in the determination of the dephasing times 
of an impulsively driven mechanical system is the avail- 
able time window of 1 ns between two successive pump 
pulses. For any coherent vibrational excitation with a 
dephasing time larger than the time-delay between two 
pump pulses, the dephasing time cannot be determined 
accurately in the time domain. Here we show how this 
limitation is circumvented by sub-harmonic resonant ex- 
citation of the phononic system with a frequency reso- 
lution achieved which would require a mechanical delay 
line of 150 m. 

The pump-probe experiments were performed at 
room temperature in reflection geometry [l2| . An 
average pump and probe power of 100 mW and 10 mW, 
respectively, is focused onto a spot of about 50 /im in 
diameter. The pump and probe wavelengths were set 
to 825 nm with a pulse length of 70 fs. The investigated 
free-standing membranes were prepared from a standard 
[100]-oriented silicon-on-insulator (SOI) wafer, which 
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FIG. 1: (Color online) Time-resolved modulation of the re- 
flected probe beam AR/Ro measured with a repetition rate 
= 800 MHz. The inset shows these oscillations after sub- 
tracting the electronic contribution (Note the change of three 
orders of magnitude in the j/-axis between main figure and 
inset.) 
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FIG. 2: (Color online) (a) Fourier spectra of the transients 
recorded with two different repetition rates. The red curve 
represents the out-of-resonance situation (/^ = 800 MHz), 
while the blue curve represents the exact resonant condition 
(/^ = 999.99MHz). Odd confined modes are labeled with 
their mode index (n = 1, 3, . . . 15). (b) Resonant behavior of 
the amplitude of fundamental mode /„ — 19 GHz (n = 1) 
when fj^ is varied as indicated. 



was produced via the smart-cut process [l5j |. The 
membranes were obtained by steps of masking, thermal 
oxidizing, and etching |13| . The single crystalline 
free-standing membranes were rectangular, roughly 
1x1 mm^ in size, and 222 nm thick. 

First we discuss the non-resonant excitation of the 
membrane, i.e. when the laser repetition rate is not 
commensurable with the fundamental acoustic mode of 
the membrane. The acoustic phonon vibrational spec- 
trum of thin free-standing membranes is characterized 
by equidistant peaks that appear as a con seq uence of 
the confinement of the acoustic modes [11, Eilll- In 
the time domain, this behavior is evidenced by the vi- 
brational superposition of a fundamental mode and os- 
cillations corresponding to odd higher harmonics [13[ . A 
transient from a pump-probe experiment taken with a 
non-resonant repetition rate of = 800 MHz is shown 
in FigUl 

The transient showing the temporal modulation of the 
reflected probe beam AR/Rg is characterized by a sharp 
rise of the signal at t = Ops due to the electronic ex- 
citation followed by an exponential decay. The features 
of interest for this work are the small oscillations that 
modulate the signal. The inset (in Fig[T]) shows these ex- 
tracted oscillations after the subtraction of the electronic 
contribution which can be modelled by a decaying multi- 
ple exponential function. The numeric fast Fourier trans- 
form (FFT) of the extracted oscillation clearly exhibits 
the narrow equidistant peaks corresponding to the longi- 
tudinal acoustic phonon modes that are confined perpen- 
dicular to the membrane's surface [Figl2^ (red curve)]. 
Due to selection rules in the generation and detection 
processes, only the odd confined modes are observed [l3| . 
The frequency of these modes is give n by /„ = 
{n = 1, 3, 5 . . .), where Vac = 8433 m/s [l7| is the longi- 
tudinal bulk acoustic velocity and d is the thickness of 



the membrane. For this membrane of d = 222 nm, the 
fundamental frequency corresponds to fi = 19 GHz [l3|, 
and modes up to n = 15 (285 GHz) are observed. 

The amplitude of the coherent oscillations barely de- 
cays within the time- window of observation of ~ 1 ns, i.e. 
the characteristic lifetime (r) of the phonon vibrations is 
much longer than 1 ns. At the arrival time of the next 
pump pulse the phonon system is still coherently excited. 
In order to impulsively excite acoustic modes construc- 
tively with subsequent pump pulses, we adjust the repe- 
tition rate of the pump laser to match the phonon oscilla- 
tions. The resonant excitation of the fundamental mode 
is carried out by tuning the repetition rate of ^ 1 GHz 
exactly to the 18*'' sub-harmonic of the membrane's fun- 
damental mode, i.e. = /i/19. In resonance, every 
19*'* tooth of the pump laser comb matches one tooth 
of the comb of the membrane's phonon mode. In the 
time domain, in this situation the phonon mode is ex- 
cited by the pump pulses once every 19*'* oscillating cycle 
of the phonon vibration, hence the modes are amplified 
resonantly. This resonance enhancement can be clearly 
observed in FiglS) The red and green curves show the 
close-to-resonance situation, whereas the blue curve rep- 
resents the exact resonance situation. The corresponding 
extracted acoustic oscillations are shown in Fig|31 The 
strong amplitude change of the generated reflection mod- 
ulations is clearly visible, as closer the pump repetition 
rate approaches the 18*'* sub-harmonic of /i. In addition 
a phase shift is observed when is varied. 

Figure [2] shows the Fourier transform of the extracted 
oscillations depicted in Fig|4^. Panel (a) of Figl2] com- 
pares the two extreme situations: the far-from-resonant 
situation (/^ = 800 MHz, red curve), and the blue curve 
corresponding to the exact resonant condition (/^ = 
999.99 MHz). Both curves show the intense and nar- 
row peaks at the modes frequency fn{n = l to 15). An 
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FIG. 3: (Color online) Typical transients for different pump 
beam repetition rates (/^) as indicated. = 800 MHz cor- 
responds to the far-from-resonance situation (also shown in 
FigH}, and = 999.99 MHz to the exact resonant condition 
with the 18*'' sub-harmonic of the membrane's fundamen- 
tal confined acoustic mode /i — 19 GHz. The transients are 
shown vertically shifted for clarity. 
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FIG. 4: (Color online) Oscillations corresponding to the mem- 
brane's confined longitudinal acoustic modes for different ex- 
citation repetition rates fa. (a) Experimantal oscillations, 
derived from Fig[3] (b) Corresponding simulated transients. 
The color legends correspond to the one shown in FigO 



increase of about a factor of 20 is observed for the n — 1 
phonon mode. The progressive increase in spectral am- 
plitude of this mode is shown for different repetition rates 
close to resonance in Fig 12)3. The width of all peaks in 
this figure is limited by the experimental time window of 
observation of 1 ns. 

The resonant excitation of the fundamental /i acoustic 
mode of the membrane is illustrated in FiglSJj. The evo- 
lution of the spectral Fourier amplitude (see Figl^Jj) cor- 
responding to this mode is plotted as function of the exci- 
tation repetition rate (bottom axis) , which is scanned 
between 997.6 MHz and 1001.6 MHz. 

The amplitude of the mode gets enhanced showing a 
Lorentzian shaped curve with its maximum exactly at 
= /i/19 = 1 GHz. The maximum amplitude of the 
membrane oscillations of AR/Rq = 1.6 x lO^'* corre- 
sponds to a change in the thickness of the membrane 
of 1.6 pm which is measured with 1 fm resolution for a 



noise limit of AR/Rq = 10~^. Furthermore, from the 
extracted oscillations in the time domain (see Fig|4^) the 
phase shift of the oscillations relative to the resonance 
condition can be derived. The evolution of this relative 
phase follows the typical phase progression of a driven 
oscillator (FiglS^). 

In order to obtain a thorough understanding of the 
membranes dynamics, we compared the experimental re- 
sults to calculations performed using a detailed model 
for the optical generation and detection processes, based 
on a macroscopic elastic and electromagnetic theory 

The large penetration depth of the exciting laser pulse 
(a-i ^ 8^m. :$> d) compared to the fast diffusion 
of the excited carriers through the thin membrane leads 
to a very fast homogenization of the carrier distribution. 
Hence, the optically generated stress can be considered 
as spatially homogeneous on our experimental time scale 
(~ 1 ns) [13| . In a similar way, since the excitation pulse 
duration (70 fs) is much shorter than the ty pica l electron- 



hole pair lifetime in silicon (10 — 100 /xs) [2l|, the tem- 
poral profile of the optically generated stress after one 
excitation can be described as a Heaviside function. The 
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FIG. 5: (Color online) (a) Phase of the extracted oscillations 
with respect to the resonant condition, and (b) normalized 
spectral amplitude of the /i acoustic confined mode as func- 
tion of the laser repetition rate . The full lines are the result 
from the simulations using the theoretical model described in 
the text. The arrows indicate the corresponding positions 
of the transients and Fourier spectra shown in the Figs[2]to|4] 
for the respective colors. The top axis indicates the frequency 
of the driven AR/Ro oscillations, i.e. 19 times 
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periodic excitation is described by additionally multiply- 
ing the response in the Fourier domain by a Dirac comb 
of frequency . Consequently the reflectivity change in- 
duced by the confined vibrations is obtained [13|. 

In the simulations, the phonon damping constant 7 
is assumed to be the same for all confined modes and is 
left free to adjust to the experimental data. The results 
for the calculated sub-harmonic resonant scan over the 
fundamental harmonic mode (n=l) are shown in FiglS] 
The simulated Fourier phase and amplitude for varying 
excitation repetition rate (/^) are plotted in panel (a) 
and (b) respectively. The best agreement with the exper- 
iment is obtained for 7 = 0.11 x lO^rad/s. This yields a 
full width at half maximum (FWHM) of the fitted curve 
of A/o = 3.2 MHz at the central frequency fo — ^ GHz, 
as depicted in Figl5][b). Since the mechanical resonance 
is at 19 times this frequency, /i — 19 fo, its FWHM is 
consequently A/i = 19 A/o = 2^/3 ^ 7 = 60.8 MHz. 

The experimental results have also been compared to 
simulations in the time domain. This is shown in Fig|4lD. 
As can be observed, the overall agreement with FigH^ 
is considerably good for the amplitude, frequency and 
the phase. For each repetition rate only the absolute 
amplitude AR/Rq for the resonant condition (blue curve) 
has been adjusted. 

In analogy to a mechanical system, the quality factor 
of the single crystalline free-standing silicon membrane 
can be defined as Q = /i/A/i = 313 for the fundamen- 
tal mode. Such a value represents a high Q-factor for a 
mechanical oscillator system, taking the frequency in the 
GHz-regime and room temperature conditions into ac- 
count. This value for Q is almost twice the one reported 
for this frequency range and temperature p^ . 

Assuming a homogeneous broadening, leads conse- 
quently to a lifetime (r = I/7) for the acoustic mode 
of r = 9 ns. This value is a lower limit for the modes life- 
time, since inhomogeneous contributions like thickness 
fluctuations are not considered in the analysis. The value 
T = 9 ns is found to be slightly larger but of the same or- 
der of magnitude as compared to those reported in earlier 
works in bulk silicon [S^, [l^l based on picoseond ultra- 
sonic experiments in this frequency range 20 GHz) . 
However, it is important to notice that the method pre- 
sented here to obtain the lifetime is conceptually more di- 
rect and accurate. In conventional picosecond ultrasonic 
experiments a metallic transducer (i.e. a thin metallic 
layer) is used to launch a hypersound pulse into the ma- 
terial, and the lifetime is indirectly determined by ana- 
lyzing the relative acoustic attenuation between the re- 
flected acoustic echoes (23 - l23 |. By the intrinsic nature 
of this method, the pulses have a more or less defined 
central frequency, but they are in general spectrally very 
broad. In addition the evaluation of the phonon lifetimes 
requires the knowledge of additional parameters such as 
the Griineisen parameter, which contributes to a large 
uncertainty in the exact numbers. By this sub-harmonic 



resonant driving method, the linewidth is directly mea- 
sured without the need of making further assumptions, 
and in contrast to Refs. 22-23] we determine the lifetime 
of a single mode. 

In conclusion, we have characterized a phononic 
system by sub-harmonic-resonant excitation of phonon 
modes using a 1 GHz femtosecond laser system. The 
accessible spectral resolution of this repetition rate 
scanning ASOPS based technique lies in the MHz range, 
which is well below the conventionally possible spectral 
resolution of measuring a time window of 1 ns. Such a 
frequency resolution is difficult to achieve with conven- 
tional pump-probe techniques, since 1 MHz resolution 
would require a mechanical scanning of a time delay 
of 1 fis corresponding to 150 m. The quality factor Q 
of the investigated membranes was determined to be 
313 at room temperature. We have established a lower 
limit for the lifetime of the membrane's fundamental 
confined acoustic mode of 9 ns. This new characteri- 
zation method enables a thorough investigation of the 
mechanical, thermal and elastic properties of thin mem- 
branes, which might be beneficial for the understanding 
and development of future nanomechanical and optome- 
chanical systems. The coherent sub-harmonic resonant 
diving might also provide an alternative, interesting 
and effective way to achieve coherent generation and 
amplification of hypersound in terahertz opto-acoustic 



devices |25l . |26| , or as well to amplify spin coherence in 



semiconductor nanostructures 
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